Evidence of kinetically arrested supercooled phases in the pervoskite oxide NdNiC>3 
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We report the time and temperature dependent response of thermopower in the non-magnetic 
perovskite oxide NdNi03. We find that on cooling below the metal- insulator transition temperature 
the system evolves into a phase separated state which consists of supercooled metallic and insu- 
lating phases. This phase separated state exhibits out of equilibrium features such as cooling rate 
dependence and time dependence. The existence of these dynamical features is attributed to the 
transformation of supercooled metallic phases to the insulating state. On cooling a small fraction of 
supercooled phases gets kinetically arrested in a glassy state and these supercooled phases remain in 
that state down to low temperature. In the heating cycle the arrested states dearrest above 150 K 
and this results in the reappearance of time dependent features. 
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I. INTRODUCTION 

First-order phase transition (FOPT) has been a topic 
of fundamental importance in the field of condensed mat- 
ter physics. The materials undergoing such a transition 
often exhibit a phase-separated state in the vicinity of 
the transition. The phase separated state in many FOPT 
compounds displays a variety of interesting phenomena 
such as colossal magnetoresistance, CMK^^ time de- 
pendent resistivity and magnetization^^ dependence 
of physical properties on rate of temperature change 
and enhanced 1// noise.— From the last decade on- 
wards, the main focus of research in FOPT studies is 
centered around understanding the origin of the phase 
separated state and its response to various perturba- 
tions such as magnetic and electric fields, pressure, x- 
ray and laser illumination, strain, temperature variation 
etc . 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 These studies reveal that coexisting 
phases have nearly equal free energies and the applica- 
tion of an external perturbation decreases the free energy 
of one phase relative to the other phase giving rise to a 
large change in the observed physical properties. Later it 
has been found that even in the absence of any apparent 
perturbations, such as mentioned above, the relative vol- 
umes of the coexisting phases can shift such that the total 
free energy decreases, which results in time dependence 
of physical properties at a fixed temperature.—^ 

The existence of a mixed phase region with nearly 
equal free energies have its origin in the phenomenon of 
metastability, in particular supercooling or superheating. 
In a first order transition where heterogeneous nucleation 
is missing, the high temperature phase can survive be- 
low the thermodynamic transition temperature Tq as a 
supercooled (SC) metastable phase.— The homogeneous 
nucleation will start when the thermal energy is of the or- 
der of the free energy barrier that separates the SC state 
from the ground state.— Each SC phase has a character- 
istic temperature T* known as the limit of metastability. 
Below this temperature the energy barrier separating the 



metastable state from the ground state vanishes and the 
SC phase is free to flip to the ground state. The SC 
phase will flip to the ground state only when their kinet- 
ics are not frozen ie their glass formation temperature 
T g is less than the limit of metastability T*M- Now if 
the glass forming temperature is higher than the limit 
of metastability, the SC phase will remain trapped down 
to low temperatures even when the free energy barrier 
vanishes. On subsequent heating the kinetics of these 
trapped phases will be restored to the metastable state 
above the temperature T > T g . In a phase separated sys- 
tem where the SC phases are not trapped in a glass-like 
state, the thermally activated stochastic switching of SC 
phases to the ground state will give rise to non-Debye 
relaxation in physical properties, 1 - So for T g < T* time 
dependent effects will be observed only in the cooling 
cycle, while for T g > T* time dependent effects will be 
observable both in heating and cooling cycle. 

Here we have used the compound NdNi03 to under- 
stand supercooling and glass formation in a polycrys- 
talline material where a distribution of SC regions co- 
exist, with each one of the different SC regions having a 
unique T*M NdNi03 belongs to the family of rare earth 
nickelates and it undergoes a temperature driven metal 
to insulator transition and a para to antiferromagnetic 
transition at 200 K.— ^ Now for the study of the dy- 
namics of a phase separated system such as NdNi03 the 
ideal probe to determine the fractional volume of the con- 
stituting phases could perhaps be a muon spin relaxation 
experiment.— But for the lack of access to such facilities 
we depend on transport properties such as resistivity and 
thermopower to probe the system. These measurements 
suffer from the effects of percolation and inter-crystallite 
contact resistance and do not allow us to do an accu- 
rate estimate of the constituent volume fractions at tem- 
peratures far below the percolation threshold^ While 
thermopower and resistivity measurements on NdNi03 
as functions of temperature has been reported earlier by 
others, in this work we focus our attention on the dy- 
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namical aspects of the M-I phase transition through time 
dependent thermopower measurements. We find that at 
temperatures far below the percolation threshold, ther- 
mopower measurements are more sensitive to the pres- 
ence of minority phases than that of resistivity. Further 
investigations suggest that on cooling the sample below 
the metal- insulator transition temperature, Tmi, while a 
majority of SC regions transform to the insulating state, 
a small fraction of SC regions get kinetically arrested in 
a glassy state. These SC regions present in the glassy 
metallic state exist down to low temperatures. 



mopower as a function of time. In the heating cycle the 
time dependent thermopower was done in a similar way: 
first take the sample to 220 K, wait for half an hour, then 
cool at 1.0 K/min to 85 K and then heat at 1.0 K/min 
to the temperature of interest, and once the temperature 
stabilized record the thermopower as a function of time. 
It was ensured that the temperature difference between 
the two ends of the sample was kept fixed at 2 K during 
all these measurements. 



II. EXPERIMENTAL DETAILS 

Polycrystalline NdNiC>3 samples in the form of 6 mm 
diameter and 2 mm thick pellets were used in the ther- 
mopower measurements. The details of sample prepa- 
ration and characterization are described elsewhere^ A 
temperature of 1000°C and oxygen pressure of 200 bar is 
required to get good quality samples. 

Below Tmi (200 K) NdNiC>3 is not in thermodynamic 
equilibrium and it relaxes with time and thus any ex- 
perimental result would depend on the procedure used 
for the measurement. We used the following procedure. 
While cooling we start with 300 K with the two sides 
of the NdNiC>3 pellet kept at a temperature difference of 
2 K, and then record the voltage (V mea s) in steps of 1 K 
after allowing the temperature to stabilize at each point 
for about 30 s. In between two temperature points the 
sample was cooled at a fixed cooling rate of 2 K/min. 
During this process also care was taken to ensure that 
the temperature difference between the two ends of the 
sample remained at 2 K.— This procedure ensures that 
the two ends of the sample smoothly vary with tempera- 
ture without any intermediate thermal oscillation. This 
condition is important since in our earlier measurements 
we have observed that thermal oscillations change the 
physical state of the system.— To remove the effect of 
any stray thermo emf in the system, we recorded the 
voltage across the sample with the ends at zero temper- 
ature difference (Vtare) in a different run under identical 
experimental conditions. This stray voltage Vtare was 
subtracted from V meas at each temperature before calcu- 
lating the absolute thermopower. After the cooling run 
was over the heating data were collected at every 1 K 
interval. The heating rate between two points was the 
same as the cooling rate used earlier. This cycle of mea- 
surement was also repeated with a different cooling and 
heating rate of 0.2 K/min. 

It was observed that the thermopower above 200 K 
does not show any time dependence, and is independent 
of the history of the measurement. So to avoid the ef- 
fect of any previous measurements all the time dependent 
measurements used the following protocol: for cooling cy- 
cle, first take the sample to 220 K, wait for half an hour, 
and then cool at 1.0 K/min to the temperature of interest 
and once the temperature is stabilized record the ther- 
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Figure 1: (Color Online) (a) Temperature versus ther- 
mopower for NdNiOs for the cooling/heating rate of 
2.0 K/min and 0.2 K/min. (b) Temperature versus resistiv- 
ity for NdNiC>3 for the cooling/heating rate of 2.0 K/min and 
0.2 K/min. Inset shows the expanded view of the resistivity 
plot, obtained while cooling, close to the temperature where 
the metallic percolation ceases (144 K). 
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III. RESULTS 
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Figure 1(a) shows the thermopower of NdNiC>3 as a 
function of temperature. The thermopower is multiple 
valued, the cooling and heating data differing signifi- 
cantly from each other and forming a large hysteresis 
loop. The sign of thermopower remains negative both in 
the metallic and insulating state which indicates that the 
charge carriers are electrons in both states. On cooling 
the sample below 300 K, the magnitude of thermopower 
decreases in magnitude till 200 K. Below this tempera- 
ture the thermopower shows a strong temperature depen- 
dence with its magnitude attaining a maximum at about 
115 K. The thermopower data indicates that NdNiC>3 un- 
dergoes a metal to insulator transition at 2 00 K which 
is consistent with the resistivity (see Figure 1(b) ) and 
earlier thermopower measurements^ Below 115 K, the 
thermopower data does not follow the band gap model 
(\S\ oc 1/T), unlike resistivity. A small jump in the 
thermopower has been observed at around 105 K, which 
is the region where the volume fraction of the metallic 
phases is small. Such behavior has been predicted from 
Bergman and Levy's theory of effective thermopower of 
mixed phase systems and has been experimentally ob- 
served in Al-Ge thin film a 24 i 27 . 

Below the M-I transition, thermopower shows a signif- 
icant d epend ence on the rate of temperature change. See 
Figure 1(a) A slower cooling rate of 0.2 K/min yields 



a higher absolute thermopower than that of the faster 
cooling rate 2.0 K/min. This rate dependence of ther- 
mopower is more readily seen below 125 K during a cool- 
ing run, while in resistivity, the cooling rate dependence 
is more significant near 145 K where the percolation of 
the metallic regions cease. In the heating cycle ther- 
mopower exhibits a large rate dependence at low tem- 
peratures and this rate dependence persists upto about 
180 K before disappearing as we increase the tempera- 
ture towards Tmi- In contrast to this we did not find 
any significant h eatin g rate dependence in the resistivity 
data(See Figure 1(b) I. The presence of cooling and heat- 
ing rate dependence in the thermopower data indicates 
that the system is not in its thermodynamic equilibrium 
both in cooling and as well as heating cycle. These ob- 
servations are corroborated by the data shown in Figure 
Hand Figured 

Figure [2] displays a subset of the time dependent ther- 
mopower data recorded at fixed temperatures in the cool- 
ing run. The data are presented as S(t)/So, where 
So = S(t — 0), so that the values are normalized for easy 
comparison. We found that below 170 K, the absolute 
value of thermopower increases with time. A maximum 
relative increase in thermopower of about 14 % was ob- 
served for a period of one hour at 152.5 K. The time 
dependence of thermopower is lower both above and be- 
low this temperature. Below 120 K, we did not find any 
detectable time dependence in the system. All the time 
dependence data were fitted to a stretched exponential 
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Figure 2: (Color Online) Time dependence of thermopower 
while cooling, at various temperatures in the range of 157.5- 
120 K, for a period of one hour 



// 


T(K) 


Si /So 


r (10 3 s) 


7 


X 2 /DOF 


R 2 


1 


140.0 


0.0426(1) 


3.12(2) 


0.729(1) 


0.34 


0.99705 


2 


145.0 


0.0812(1) 


3.09(2) 


0.782(1) 


1.3 


0.99837 


3 


147.5 


0.0905(2) 


2.53(1) 


0.896(1) 


2.3 


0.99853 


4 


150.0 


0.118(1) 


2.34(1) 


0.826(1) 


3.2 


0.99861 


5 


152.5 


0.231(1) 


3.62(2) 


0.776(1) 


5.3 


0.99905 


6 


155.0 


0.0673(2) 


2.58(1) 


0.834(2) 


7.2 xlO 3 


0.99687 



Table I: Fit parameters for a subset of time dependence data 
shown in Figure [2] The degrees of freedom of the fits DOF ~ 
1000. The x 2 / D OF for 155 K is too large, indicating that we 
have underestimated the error in thermopower in this case. 
Anyway, we note that, the R 2 values are consistently good 
and indicate reasonably good fits. 



form 



S(t) = S Q + S 1 (l-e-(*) 7 ) 



(1) 



where Sq, S%, t and 7 are fit parameters with S(t), So and 
Si having negative value. The fits are quite good with 
R 2 value greater than 0.995 in most cases. See table 1. 

Figure [3] shows a part of time dependent thermopower 
data recorded at fixed temperatures in the heating runs. 
We did not find any significant time dependence up to 
150 K. Above 150 K, the magnitude of thermopower in- 
creases with time. The maximum relative increase in 
thermopower in one hour is around 1% which is much 
smaller than that seen in the cooling runs. Above T MI , 
the thermopower data is stable and no time dependence 
or cooling/heating rate dependence was observed . 
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Figure 3: (Color Online) Time dependence of thermopower 
while heating from 85 K, at various temperatures in the range 
of 150 -200 K, for a period of one hour 



IV. DISCUSSION 

The presence of out of equilibrium features in ther- 
mopower is similar to that observed in our resistivity 
measurements reported earlier^ Based on those obser- 
vations we had proposed a model to understand the be- 
havior of the phase separated regime. The model is as 
follows: The high temperature metallic phase survives 
below the metal to insulator transition in its supercooled 
metastable state. The supercooled metallic phase is not 
a single entity, but instead is made up of various regions 
which can make a transition or switch from the metal- 
lic to the insulating state independent of the neighboring 
regions. These regions are referred to as switchable re- 
gions (SR). Each SR in its metallic state has to cross an 
energy barrier U to go to its insulating state and the bar- 
rier strength was determined to be U — cV(T — T*) 1 / 4 , 
where V is the volume of the SR, T* is its temperature 
of metastability, and c is a constant. An SR will switch 
to the insulating state if (a) it attains its temperature of 
metastability or (b) it gains enough energy from thermal 
fluctuations so that it can cross the energy barrier U. 
For T* <T < T M i and T g < T, supercooled phases will 
switch over to the insulating state stochastically which 
will cause time dependence in the physical properties. 
No time dependence in resistivity was observed in the 
heating cycled On the basis of this it was claimed that 
all the SR's switch to the insulating state on cooling im- 
plying that T* > T g for all of them. 

In thermopower measurements, during cooling, we got 
the maximum time dependence at around 152.5 K while 
in the case of resistivity the maximum time dependence 
occurs around 145 K. The maximum magnitude of change 
in thermopower in one hour (AS) is ~ 14 %, which is 
much smaller than that of resistivity (Ap ss 200%). 17 



We notice that the change in thermopower in one hour 
AS is much closer to the change in insulating volume 
fraction (AV « 8%).— The time dependence of ther- 
mopower in the cooling cycle can be understood on the 
basis of our model. As we have discussed in the previ- 
ous paragraph, below Tmi, the SR's present in the SC 
state stochastically switch to the insulating state which 
enhances the volume fraction of insulating phase with 
time. This increase in insulating volume fraction results 
in a higher magnitude of thermopower. Thermopower 
being a transport property, the change in thermopower 
is also affected by percolation effects and thus the be- 
havior of time dependence of thermopower and the time 
dependence of insulating volume fraction will not have a 
one-to-one monotonic correspondence. 

Below 120 K no detectable time dependence was ob- 
served in the thermopower, but cooling rates of 2 K/min 
and 0.2 K/min lead to slightly different low temperature 
states with differing thermopower. This suggests that 
a tiny fraction of SR's must remain in the supercooled 
state even below 120 K. Since we did not find any time 
dependence below 120 K, we argue that the kinetics of 
these SR's are arrested and and they end up in a glass- 
like state. Now, the question arises why there is no rate 
dependence in resistivity below 120 K neither in cool- 
ing runs nor in the subsequent heating runs. To answer 
this question we suggest that the small fraction of SR's 
trapped in the glassy state are not detectable in a re- 
sistivity measurement. This idea is in consonance with 
the fact that thermopower measurements show a small 
time dependence of about 1% over a period of an hour 
in heating runs (Figure [3| while no detectable time de- 
pendence was observed in the heating runs of resistivity 
measurements.— 

The advantage of thermopower over resistivity in de- 
tecting tiny amounts of SC metallic regions (SR's) dis- 
persed in an insulating matrix can perhaps be understood 
as follows. When the fraction of SC metallic regions is 
small they will be completely surrounded by insulating 
regions. Now the net contribution of a tiny metallic re- 
gion embedded in an insulating matrix on the total resis- 
tivity of the material will be determined by the resistance 
of the metallic region plus its inter-crystallite resistances 
in the current path. Now we conjecture, based on our 
experience on the contact resistance of polycrystalline 
NdNiOs, that the inter-crystallite resistance between a 
metallic and an insulating crystallite will be very high 
compared to the resistance of the metallic crystallite. 
This implies that the tiny metallic region will be essen- 
tially shut out of the current path and will not be able 
to affect the resistivity of the material significantly. The 
boundary effects are much less significant in the case of 
thermopower compared to resistivity because the thermal 
conductivities of the metallic and insulating regions are 
comparable and this makes the thermopower sensitive to 
the embedded metallic regions . 28 ' 29 

In a heating run, the SR's present in a glassy state will 
remain arrested in that state till their respective T g 's 
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are attained. In cooling runs we found time dependence 
down to 120 K while in heating runs we did not find 
any time dependence until 150 K which suggests that 
the SR's which are arrested in the glassy state have their 
T g 's above 150 K. On heating the system above 150 K, 
the kinetics of the SR's gets dearrested and it can cross 
the free energy barrier U and flip to the insulating state 
with a probability of flipping oc exp(—U /ksT). In the 
time dependence data presented in Figure [3] the flipping 
of the SR's to the insulating state raises the insulating 
volume fraction with time and this enhances the magni- 
tude of thermopower with time. The fact that we get 
noticeable time dependence all the way upto Tmi sug- 
gests that the T g 's are distributed quite broadly upto 
Tmi ■ We might now say that each SR has a unique tem- 
perature of metastability T* and a unique temperature 
of kinetic arrest T g . In a cooling run, the majority of the 
SR's present in the supercooled state switch to the insu- 
lating state by the time their T* is attained and only a 
small fraction of SR's get arrested in the glassy state and 
this suggests that the majority of the SR's have their T* 
greater than T g . 

V. CONCLUSION 

Our experimental results and analyses suggest the fol- 
lowing picture for the NdNi03 system. Below Tmi the 



system consists of SC metallic and insulating regions. 
Each of the SC metallic region has a unique tempera- 
ture of metastability T* and a unique temperature of 
kinetic arrest T g . For the majority of the SC regions T* 
is greater than T g . Thus when the sample is cooled be- 
low Tmi, while the majority of the SC regions switch to 
insulating state when their temperature of metastabil- 
ity is reached, a small fraction of the SC regions which 
have T g > T* get arrested in a glassy state and remain 
in that state down to low temperature. Such SC regions 
have their T g greater than about 150 K. In a heating run, 
these SC metallic regions remain trapped in their glassy 
metallic state till their temperature of kinetic dearrest 
(T g ) is reached. Above this temperature, these regions 
switch over to the insulating state stochastically. 
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